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ABSTRACT 

Ultra-compact dwarf galaxies (UCDs) are stellar systems with masses of around 10^ 
to 10^ M0 and half mass radii of 10-100 pc. They have some properties in common 
with massive globular clusters, however dynamical mass estimates have shown that 
UCDs have mass-to-light ratios which are on average about twice as large than those of 
globular clusters at comparable metallicity, and tend to be larger than what one would 
expect for old stellar systems composed out of stars with standard mass functions. 

One possible explanation for elevated high mass-to-light ratios in UCDs is the 
existence of a substantial amount of dark matter, which could have ended up in UCDs 
if they are the remnant nuclei of tidally stripped dwarf galaxies, and dark matter 
was dragged into these nuclei prior to tidal stripping through e.g. adiabatic gas infall. 
Tidal stripping of dwarf galaxies has also been suggested as the origin of several mas- 
sive globular clusters like Omega Cen, in which case one should expect that globular 
clusters also form with substantial amounts of dark matter in them. 

In this paper, we present collisional N-body simulations which study the co- 
evolution of a system composed out of stars and dark matter. We find that the dark 
matter gets removed from the central regions of such systems due to dynamical fric- 
tion and mass segregation of stars. The friction timescale is significantly shorter than 
a Hubble time for typical globular clusters, while most UCDs have friction times much 
longer than a Hubble time. Therefore, a significant dark matter fraction remains within 
the half-mass radius of present-day UCDs, making dark matter a viable explanation 
for the elevated M/L ratios of UCDs. If at least some globular clusters formed in a 
way similar to UCDs, we predict a substantial amount of dark matter in their outer 
parts. 
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1 INTRODUCTION 



Ultra-compact dwarf galaxies (UCDs) were discovered in 
the late 1990s in spectroscopic surveys of the Forn ax galaxy 
cluster l|Hilker et al.lll999l : lDrrnkwater et al.ll2000l) and have 
since th en been found in other ne arby galaxy clu sters as well 
(,Haseganet all 120051: iMieske et al.. ,2005 . .20071 : Ijones et all 



20071 ). They are 



l2006l : jPirth et all l2007l : iReikuba et al 
bright (-11 < Mv < -13.5) and compact (7 < < 100 pc) 
stellar systems which have ages of at least several Gyr and 
possi bly up to 10 Gyr (|Mieske et al.|[2006l : lEvstigneeva et al.l 
l2007l ). The masses and sizes of UCDs are larger than those 
of Galactic globular clusters, but similar to those of nuclei 
in dwarf elliptical galaxies (Drinkwater et al. 2003, Bekki et 
al. 2003). 

One of the most remarkable properties of UCDs is 
that their dynamical mass-to-light ratios are on average 



about twice as large than those of globular clusters of 
comparable metallicity, and also tend to be larger than 
what one would expect based on simple stellar evolu- 
tion models that a ssume a stand ar d stellar initial mass 
function, like e.g. Kroupal ( 2001 ) jHasegan et al. 20051 : 



iDabringhausen. Hilker fc KroupalHoOSi : iMieske et al.|[20ci8l ) 
If not due to a failure of stellar evolution models, this points 
either to unusual stellar mass functions (IMieske fc Kroupal 
l2008l : [Dabringhausen. Baumgardt fc Kroupall200i ^ or possi- 
bly to the presence of a significant amount of dark matter in 
UCDs. We note that most methods used to determine M/L 
ratios for UCDs rely on the assumptions that mass follows 
light and isotropic velocity dispersions. How well these as- 
sumptions are fulfilled is currently not known. Also due to 
the large distances of UCDs, only integrated velocity disper- 
sions can be obtained, which in most cases are intermediate 
between the central and the global velocity dispersions. In 
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both GCs and UCDs are formed with the same non-zero 
dark-to-stellar-mass-fraction. We then investigate how the 
dynamical co-evolution of dark matter and stars changes 
the observed dark matter fraction as a function of time. 
We assess whether the observed rise of M/L ratios from the 
regime of GCs to that of UCDs can be explained by our 
working hypothesis and the subsequent dynamical evolution. 
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Figure 1. The mean mass density within the half-mass radius of 
the joint sample of GCs and UCDs from Fig.[5]is plotted vs. their 
relaxation time llMieske et al .120081) . The dotted line indicates the 
approximate central (r < lOpc) da rk matter densities expected 
for cuspy dwarf galaxy CDM halos llCilmore et al.ir2007ft . 



order to determine the mass-to-light ratio, the mass mod- 
eling has to take the density profiles of the UCDs as well 
as the effects of seei ng and a fini t e slit size into account, as 
done for example in lHilker et al.l l|2007l ). 

Several formation scenarios have been discussed for 
UCDs, like e.g. that UCDs are simpl y massive globular 
clusters and form in t h e same way llHilker et al.l 1 19991 : 
lEvstigneeva et all l2007l : iForbes et all 120081 ). that they 
are the nuclei of tid ally stripped, originally muc h more 
extended galaxies 
Bekki et al. . 2003. ; 



('Bek ki. Couch fc Dr inkwatei j l200ll: 



^ Thomas. Drinkwater fc Evstigneeval 

20081 : ICoerdt et al.ll2008l), or t hat they are merged glob ^ 



ular clusters iOh fc LinI |200GI : iFellhauer fc Kroupal l2002h . 
ICoerdt et al.l 1 20081 ) have shown that tunneling of dark 
matter to the central region of a disk galaxy, due to 
gas-infall, can significantly increase the M/L ratios in 
the rmclear region, and hence may explain the elevated 
M/L ratios of UCDs, provided that UCDs formed by tidal 
stripping. Indeed, it has been suggested that also GCs 
may have originated as centers of individual primordial 
dark matter halos (e.g. Carraro fc Lia 2000, Lee et al. 
2007, Bekki et al. 2007). If dark mat ter tunneling is an 
efficient mechanism (iGoerdt et ahlliooil ). one may therefore 
expect both UCDs and GCs to be formed with a significant 
fraction of dark matter. It is important to note that such an 
increase of dark matter density by some kind of tunneling 
mechanism is necessary to explain a significant amount of 
dark matter in UCDs or GCs, since their present-day stellar 
(and hence implied dark matter) densities are up to 2-3 
orders of magnitude higher than ex pected for cuspy dar k 
matter halos of dwarf galaxy mass (|Gilmore et al.l l2007^ . 
This is shown in Fig. [1] 

In this paper, we start from the working hypothesis that 



2 THE MODELS 

In our simulations, we assume that stars and dark mat- 
ter particles follow th e same density d istribution initially, 
which was given by a IPlummeil l|l91ll ) model. Determina- 
tions of mass-to-light ratios of globular clusters or UCDs 
rely mainly on stars located inside the half-mass radius 
(Hilkcr et alll2007l: iM ieske et al."2008l) or even closer within 
( McLaug hlin fc van dcr Marcl 2005;) . Tidal effects are there- 
fore not likely to have a strong infiuence on determined mass- 
to-light ratios. Also, due to their high mass and correspond- 
ing large dissolution times, tidal interactions probably play 
only a minor role for UCDs. In our simulations we therefore 

neglect the influence of an external tidal field. 

W e assume that the stars initially follow a iKroupal 
(|200ll ) mass function with lower and upper mass limits of 0.1 
and 100 M©. Stellar evolution changes the mass function of 
stars, however most of this change happens within the first 
10^ yr s (see e.g. the grid of models by|Baumgardt fc Makincl 
(|2003h ). i.e. on a timescale short against the lifetime of glob- 
ular clusters or UCDs. We therefore also neglect stellar evo- 
lution and immediately transform stars to the assumed age 
of GCs and UCDs, T=12 Gyrs. For the transformation, we 
assume that stars with mass larger than 25 Mq form black 
holes and assume that the black hole mass is 10% of the 
mass of the initial star. This way, black hole masses in our 
models are com patible w i th ob served masses for stellar mass 
black holes (e.g. lCasared (|2006l )). Stars with masses between 
8 M0 and 25 Mq are assumed to form neutron stars with a 
mass of m = 1.3 Mq l|Thorsett fc Chakrabartvll 19991 ). Stars 
between 0.8 M0 and 8 Mq are assumed to form white dwarfs 
due to stellar e volution. The masses of the white dwarfs are 
obtained from iKalirai et al.l (|2008l ). who found, based on 
observations of white dwarfs in star clusters, the following 
relation between the initial and final mass of white dwarfs: 
rriwd ~ 0.109m -I- O.394M0. Stars less massive than 0.8 Mq 
are still on the main sequence and we assume that they have 
not yet lost any mass. The following table summarises our 
initial-to-final mass relation: 



m, m < O.8M0 

0.109 -I- 0.394, O.SMq < m < 8Mq 
1.35, 8M0 < m < 25Mq 
0.1m, 25M0 < m 



(1) 



Neutron stars and probably also black holes receive 
kicks at the time of their birth due to asymmetric super- 
nova ex plosions. The size of t hese kicks is a few hundred 
km/sec l|Lvne fc Lorimeilll994l ). which is large enough that 
most will be lost from globular clusters or UCDs (jPfahl et al.l 
I2OO2I ). We therefore assume only a small neutron star and 
black hole retention fraction of 20% in our simulations. With 
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these assumptions, the mean mass of stars in our models is 
0.344 Mq. 

The dark matter is also modeled as point mass particles. 
In our reference simulation we assume a mass of 0.03 for 
the dark matter particles, but we also make simulations with 
masses of 0.15 M© and 0.01 M© to study the influence of 
the adopted particle mass on our results. Theoretically, the 
dynamical friction of stars should not depend on the mass of 
the dark matter particles as long as th e mass ratio between 
stars and dark matter is high enough l|Binnev fc Tremaind 
1 19871 ). Also, with the adopted masses, the self-interaction 
of the dark matter particles is still unimportant over the 
timescales studied here. We will further investigate the in- 
fluence of the mass of the dark matter particles in sec. 13.41 

All runs ar e performed with the coUisional TV-body 
code NBODY4 (lAarsethl ll999l ) on the GRAPE6 computers 
l|Makino et al., 20031 of Bon n University and the results are 
expressed in A^-body units l|Heggie fc Hutll2002l ). in which 
the constant of gravity and total cluster mass are equal to 1 
and the total potential energy is equal to -0.5. Table 1 gives 
an overview of the runs performed. 



3 RESULTS 

3.1 Mass segregation timescale: analytical 
estimate 

In the following we derive an analytical estimate of the mass 
segregation time scale in compact stellar systems. 

Massive stars will segregate against dark matter parti- 
cles and lighter stars as a result of dynamical friction and 
energy equipartition. Since the masses of stars are much 
higher than the mass of the dark matt er particles, the fric- 
tional drag on the stars is given by fsee lBinnev fc Tremaind 
l| 19871 ) Eq.7-18): 



dv 
IE 



-47rln AG^p(r)m 



erf (a) r=e 



(2) 



where p(r) is the background density of dark matter and 
stars, m the mass of an inspiraling star. In A the Coulomb 
logarithm, and X = «/(\/2cr) is the ratio between the veloc- 
ity of a star and the (ID) stellar velocity dispersion o. If we 
assume z; « cr, it foll ows that X = l/y/2. Settin g InA = 12 
for globular clusters (|Binnev fc Tremaind (|l987l ) Tab. 7-1), 
the eq. [2]can be rewritten as: 



(3) 



^ = _29.96^^4^. 
dt 

The resulting energy change is ^ = ^(^^■'^^) = ''^''^TE' ^'-'^ 
a distribution of stars in virial equilibrium, an energy change 
dE corresponds to a change in potential energy md^ — 2dE. 
Hence 

^ = -59.93^!^ . (4) 
dt V ^ ' 

For a Plummer model, the density p, circular velocity Vc and 
specific potential $ at point r axe given by: 



pir) = 
<I>(r) = 



a + r 



4tt 

(„2+ ^2)1/2 



(5) 



Vc{r) 



VGMTotr^ 



,3/4 



where Mrot is the total cluster mass and a is the scale ra- 
dius of the Plummer model. With these equations, the above 
relation can be rewritten as: 



? = -14.31 

dt 



(6) 



For an order of magnitude estimate of the inspiral time scale, 
one can approximate (a^-f r^)^''* ~ a^^^. One can then 
solve the above relation and obtains as time which a star 
starting at radius Ro needs to reach the centre: 



TpHc = 0.023 Rl 



(7) 



For a Plummer model, a = 0.766-R£f , so the inspiral time for 
stars near the half-mass radius is given by: 



Tmc = 0.035 



For a globular cluster or UCD we thus obtain: 



Tmc = 5.86 



108 Mf 



1/2 



Rh 
5pc 



3/2 



m 



Gyr 



(8) 



(9) 



The resulting dynamical friction time agrees to within 20% 
with what Binney fc Tremaine found for the inspiral time 
of an isothermal sphere (their eq. 7-26). Eq. [5] predicts a 
dynamical friction time scale of 4-5 Gyr for a typical GC 
(3x10^Mq, rh = 3pc), and about 400 Gyr for a typical 
UCD (IO'^Mq, Th = 20pc). That is, after a Hubble time 
most of the dark matter in globular clusters should have 
been pushed out of the centre, while in UCDs the inspiral of 
stars should be far from complete and a significant fraction 
of DM should still reside in their centres, leading to high 
mass-to-light ratios. 



3.2 A-body results 

Fig. [2] shows the evolution of Lagrangian (lag.) radii, i.e. 
radii which contain a certain fraction of the total mass, of 
stars and dark matter particles in our first simulation, which 
had a dark matter particle mass of m = 0.03 Mq and an 
equal amount of mass in stars and dark matter. The effect 
of dynamical friction and mass segregation is clearly visible 
since the lag. radii of dark matter particles increase with 
time while those of the stars shrink. In A-body units, the 
total cluster mass and constant of gravity are both unity. 
With a mean stellar mass of m = 6.0 • 10~*, we predict 
a dynamical friction timescale of Tpric ~ 391 in A-body 
units according to eq. |S] It can be seen that by this time the 
cluster center is indeed nearly free of dark matter: inside the 
lag. radius of 10% of the stars, only 1% of the dark matter 
particles are located. At the end of the simulation, the 10% 
lag. radius of the dark matter particles is almost equal to 
the half-mass radius of the stars, i.e. only 20% of the cluster 
mass is still made up of dark matter inside the (visible) 
half-mass radius of the cluster. The core of the cluster is 
even stronger depleted and is nearly free of dark matter by 
the end of the simulations. It would therefore be difficult 
to detect the remaining dark matter by its effect on stellar 
velocities if mainly stars from the cluster center or inside 
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Figure 2. Evolution of Lagrangian radii, i.e. radii which contain 
a certain fraction of the total mass, of stars (red solid lines) and 
dark matter particles (blue dashed lines) in the first run from 
Table 1. The left panel depicts the evolution as a function of A'^- 
body time, the right panel as a function of current relaxation time, 
where the relaxation time is calculated based on the distribution 
of stars only. Inside the half-mass radius of the cluster, less than 
20% of the total mass is made up out of dark matter after the 
clusters are ten apparent relaxation times old. 



the half-mass radius are used to determine the line-of-sight 
velocity dispersion. 

The right panel of Fig. [2] depicts the evolution of lag. 
radii as a function of the ratio of cluster age to the actual 
relaxation time. In order to allow for a better comparison 
with observations, the relaxatio n time is calcu lated from the 
stellar component according to (ISpitzeilll987l '): 



Trh = 0.138 



Gm, In 7A''* 



(10) 



where M, is the total stellar mass of the cluster, Rh* 
the half-mass radius of the stellar distribution, m, and 
A^* are the mass and number of stars and 7 a constant 
in the Coulomb logari thm which is taken to be 7 — 0.11 
l|Giersz fc Heggiel ll994h . Eq. [TD]would be the relaxation time 
inferred by an observer who can only determine the stellar 
distribution and does not know about the dark matter. It 
has the same dependence on cluster mass and radius as the 
friction timescale and can therefore also be used to judge 
the dynamical state of a cluster. 

The right panel of Fig. [2] shows that once a cluster is 
two to three apparent relajcation times old, the centre is free 
of dark matter and by the time the cluster has become ten 
relaxation times old, there is little dark matter left inside 
the half-mass radius. Since most globular clusters have re- 
laxation times of only a few Gyr, their mass-to-light ratios 
should be within 20% of those of pure stellar populations if 
mainly stars inside the clusters half-mass radius are used to 
determine the velocity dispersion. Since this is within the un- 
certainty of measured mass-to-light ratios and current stellar 
population models, such a small dark matter cannot be de- 
tected kinematically in globular clusters. UCDs on the other 
hand have relaxation times significantly larger than a Hub- 
ble time and should therefore still have large mass-to-light 
ratios if they formed as a mix of dark matter and stars. 

This is confirmed by the upper panel of Fig. O which 
depicts the dark matter fraction inside the cluster core (as- 
sumed to be the region inside the 5% lag. radius of the stars) 
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Figure 3. Dark matter fraction (upper panel) and average mass 
of stars (bottom panel) as a function of time. Core values are 
shown by solid lines, values inside the half-mass radius by dashed 
lines. The core radius is assumed to be equal to the 5% lag. radius 
of the stellar distribution. The dark matter is depleted from the 
centre within 1 to 2 friction times. At the same time, heavy mass 
stars segregate against the light mass stars and the average mass 
of stars increases in the centre. Clusters where the centers are 
depleted of dark matter should therefore also be mass segregated 
in their centre. 



and inside the half-mass radius of the cluster stars. It can 
be seen that after about 1.5 to 2 friction times, the cluster 
core is almost completely free of dark matter. Within the 
half-mass radius, only about 30% of the initial dark matter 
amount remains after this time. The lower panel of Fig. [3] 
depicts the evolution of the average mass of stars in the core 
and inside the half-mass radius. At both radii, the average 
mass of stars is increasing since, while stars segregate against 
the dark matter particles, heavy-mass stars also segregate 
against the lighter ones. After about 2 dynamical friction 
timescales, the mass of stars has reached a near constant 
value of about 0.65 M©, which is significantly higher than 
the average mass of stars. Clusters which have expelled dark 
matter out of their centres should therefore also be mass seg- 
regated. 

3.3 Comparison with observations 

Fig. 2] shows the effect which the decreasing dark matter 
fraction in the center has on the projected velocity disper- 
sion of stars. In order to determine this effect, we first cal- 
culate the velocity dispersion profile (Jobsir) of bright stars 
with masses in the range 0.6 < m < 0.9 M© as a func- 
tion of projected radius. We restrict ourselves to this mass 
range since in a globular cluster or UCD, these would be 
the stars which dominate the cluster light. After determing 
the velocity dispersion profile of bright stars, we calculate 
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the expected velocity dispersion profile based on the stellar 
density distribution according to (Binney & Tremaine 1987, 
eq. 4-54): 



1 

p{r) 



Pir", 



dr 



dr' 



(11) 



where p(r) is the (3D) density distribution of bright stars 
and <I>(r) is the potential coming from the stars alone. Eq. 
111! assumes a spherical cluster potential and an isotropic 
velocity dispersion of stars. After projecting a{r) we can 
calculate the correction factor / needed so that the pre- 
dicted velocity dispersion matches the true velocity disper- 
sion of the clusters in the A^'-body simulations, i.e. /(r) = 
o'Obs(r)/apredir). The resulting correction factor is plotted 
in Fig.|4]for the first run from Table 1. Initially, dark matter 
and stars follow the same density distribution, so the veloc- 
ity dispersion is a factor f{r) = ^ {Md m + M, ) /M, = 1.41 
higher than predicted by eg. 1111 As the cluster evolves, dark 
matter is removed from the center, so the velocities of stars 
in the center are determined more and more by the stars 
alone and / approaches unity. After 3 relaxation times, the 
central velocity dispersion is only 10% higher than what 
one would expect based on the stars alone and after 10 re- 
laxation times the difference is less than 1%. Most globular 
clusters should therefore have central M/L ratios which are 
close to those predicted by stellar population models. Be- 
yond 10 half-mass radii, / remains close to the initial value 
even after 10 relaxation time s. As long as dark matte r is 
not removed by tidal effects (iMashchenko &: Sillsll200"5l ). it 
should therefore be detectable in globular clusters through 
the observation of stellar velocities in the outer cluster parts. 

We finally discuss the influence of the dark matter on 
the global mass-to-light ratios. In order to compare our sim- 
ulations with observed clusters, we again calculate true and 
expected velocity dispersions of stars with masses in the 
range 0.6 < m < 0.9 Mq. Since mass-to-light ratios of UCDs 
are determined from stellar velocities covering a significant 
fractio n of the cluster area (see e.g. discussion in lHilker et al.l 
and measured mass-to-light ratios of globular clus- 
ters are based mainly on stars in the inner cluster parts 
l|McLaughlin fc van der Marell |2005| ) , we determine global 
velocity dispersions in the simulations for all stars located 
inside the projected half-light radius. The resulting mass-to- 
light ratios of our model clusters are then given by 



M/L = f M/L\^ 



(12) 



where M / L\, is the mass-to-light ratio which a pure stellar 
population would have and / is again / = oobs/'^Pred- 

Fig. [5] depicts the evolution of M/L with cluster age for 
runs 1 and 2 of Table 1, and compares it w i th obs erved M/L 
ratios of UCDs and GCs from lMieske et all (|2008| ). Note that 
the literature M/L estimates are normalised to the same 
(solar) metallicity, to allow direct intercomparison. Time is 
again expressed in terms of age divided by the relaxation 
time as determined from the stars alone. 

The observed normalised M/L ratios show a clear trend 
in the sense that dynamically more evolved systems have 
on average lower M/L values. The mass-to- light ratios in 
our simulations also decrease as the dynamic age increases, 
since, as the dark matter is depleted from the cluster cen- 
ters, the velocity dispersion is determined more and more by 



1.4 



b 1 2 

II ^-"^ 




lJ_ 



Li_ 



0.1 1 10 

Radius [Projected half-light radii] 

Figure 4. Evolution of the ratio / of observed velocity dispersion 
to predicted one based on the stellar distribution alone for 4 dif- 
ferent ratios of cluster age to apparent relaxation time calculated 
according to eq. 1101 Initially, dark matter and stars follow the 
same distribution and there is an equal amount of mass in dark 
matter and stars, so crobs = v^fpred independent of radius. As 
the dark matter is expelled from the centre, stellar velocities in 
the inner parts are increasingly determined by the stars alone 
and / drops towards unity in the cluster centre. Globular clusters 
should therefore have central mass-to-light ratios in agreement 
with stellar population models. 



the stars alone, so M/L approaches A//L|,. Depending on 
whether a stellar M/Lj* of 2.5 or 2.0 is assumed, a run with 
a primordial dark matter content equal to or twice as high as 
the stellar mass provides an acceptable fit to the data, mak- 
ing dark matter a viable alternative to explain the elevated 
mass-to-light ratios of UCDs. We note that if the observed 
decrease of M/L with dynamical age is due to the dynamical 
depletion of non-luminous particles from the cluster centers, 
the dark matter particles have to be of lower mass than the 
stars, ruling out e.g. a central concentration of black holes 
in UCDs as the explanation for their high M/L ratios. 

The stellar mass-to-light ratios we have to assume 
in order to flt the data for GCs are marginally lower 
than predicted by simple stellar population models for 12 
Gyr old, solar-met a llicity star clusters. For example, the 
iBruzual fc CharlotI (|2003l 'l models predict an M/L of 2.5 
for a 8-9 Gyr old stellar population. The difference to a 
12 Gyr old population (M/L=3.5) is still within individ- 
ual error bars for the literature estimates, and there is also 
some uncertainty in the underlying stellar mass functions 
and stellar population models. Nevertheless, we note that 
the slightly too low M/L ratios may also be interpreted as 
signs of preferentia l loss of lo w-mass st ars in the Gala ctic 
tidal field (.Baumgardt fc Makino 200 iKruiissenI 120081 '). If 
this was the case, then less dark matter would be needed 
to explain the elevated mass-to-light ratios of UCDs, imply- 
ing a DM mass of ~50-80% of the stellar mass. However, it 
is unclear whether the actual dissolution times of the GCs 
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Figure 5. Mass-to-light ratios of UCDs (circles) and globular 
clusters (triangles) as a function of their age divided by their 
relaxation time. There is a clear trend towards lower M/L values 
for dynamically more evolved systems. The red solid and blue 
dashed curves show predicted M/L values for two of our runs 
calculated from the velocity dispersion of bright stars inside the 
clusters half-mass radius, assuming stellar mass-to-light ratios of 
M/L\t = 2.0 and M/L|* = 2.5 for the two runs. It can be seen 
that the resulting theoretical curves provide a good fit to the 
combined globular cluster/UCD sample. 



with available M/L measurements ar e short enough to 
experienced significant evaporation (iMieske et al.l 12008 



to have 
A 

case-by-case analysis for Galactic GCs will be necessary to 
assess this, based o n measured absolu te proper motions and 
orbital parameters (|Allen et al.ll200^ ). 



3.4 Scaling issues 

We finally discuss a possible biasing of our results due to the 
finite mass of the dark matter particles. Fig. |6] depicts the 
evolution of Lagrangian radii in simulations with different 
dark matter particle masses. All simulations had an equal 
amount of mass in stars and dark matter and the mass of 
the dark matter particles was set to be m = 0.1 Mq (blue 
lines), m — 0.03 Mq (red lines) and m — 0.01 Mq (green 
lines) . Since the mass of heavy stars which drive the inspiral 
is in all cases much higher than the mass of the dark matter 
particles, eq. |8] should still apply for the inspiral timescale. 
In all three simulations, the mass of individual stars if ex- 
pressed in A'^-body units was the same, so according to eq. 
[S] the inspiral timescale of the stars should be the same in 
the three simulations. 

It can be seen that the inspiral of stars and the ejection 
of dark matter particles happens in all three clusters in a 
very similar way. The agreement is especially good between 
the two simulations with the lightest dark matter particles. 
We therefore conclude that the adopted mass m = 0.03 M0 
for the dark matter particle does not influence the results 
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Figure 6. Evolution of Lagrangian radii of stars (left panel) and 
dark matter particles (right panel) in simulations which assume 
dark matter particle masses of m = 0.1 Mq (blue lines), m = 0.03 
M0 (red lines, the default mass for the simulations presented in 
Figs. [2] to [Sjl and m = 0.015 Mq (green lines). The agreement 
between the different curves, especially for light dark matter par- 
ticles is very good, showing that the adopted mass m = 0.03 Mq 
of the dark matter particle should not influence our results. 



presented in Figs. [5] to [S] Our simulations should therefore 
give a correct picture of the dynamical ejection of dark mat- 
ter from the centers of globular clusters and UCDs. 



4 CONCLUSIONS 

We have performed collisional A-body simulations of the 
evolution of compact systems composed out of a mix of stars 
and dark matter particles. Our simulations show that dark 
matter is depleted from the centers of these systems due to 
dynamical friction and energy equipartition between stars 
and dark matter particles. The inspiral time of stars is short 
enough that only 20% of the original dark matter would 
remain within the half-mass radius in typical globular clus- 
ters. If mainly stars from the inner cluster parts are used 
to determine mass-to-light ratios, the resulting increase in 
the mass-to-light ratio is within the errors with which mass- 
to-light ratios are typically determined for globular clusters 
and would therefore be difficult to detect. 

If not tidally stripped, dark matter should also reside in 
the outer par ts of globular clusters. For a number of glob- 
ular clusters, IScarpa et al.l (|2007. ) have indeed reported a 
fiattening of the velocity disperion in the outer cluster parts. 
This could however be due to a number of reasons like con- 
tamination of the sample by background stars or the tidal 
interaction of a star cluster with the gravitational field of 
the M fiky Way (jPrukier et al.lll998l : [Capuzzo Dolcetta et all 
I2OO5I ). Detailed simulations would be necessary to exclude 
these possibilities and confirm that the observed fiattening 
is due to a dark matter halo. 

UCDs on the other hand have inspiral times signifi- 
cantly longer than a Hubble time and therefore still con- 
tain most of the dark matter in their centers. Dark matter 
therefore seems a viable explanation for the elevated M/L 
ratios of UCDs, provided that UCDs originate from the cen- 
ters of dark matter halos and have seen their dark matter 
content being increased by dark matter fun neling, through 
e.g. adiabatic gas infaU l|Goerdt et al.ll2008l ). 

A prediction of our simulations, which can in principle 
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be tested by observations, is that globular clusters which 
have expelled the dark matter from their centers should 
also be mass segregated. Non-mass segregated clusters with 
velocity dispersions and mass-to-light ratios in agreement 
with simple stellar population models, would therefore have 
formed without significant amounts of dark matter in their 
centers. 

Also, if dark matter existed in a globular cluster at the 
time of its formation, it should still reside in its outer parts, 
especially if tidal stripping due to external tidal forces from 
the host galaxy l|Mashchenko fc Sillj 120051 ) and relaxation 
driven internal mass loss was not important for the clus- 
ter evolution. In this case, the measured mass-to-light ra- 
tio should increase towards the outer cluster parts, which 
can in principle be detected with dedicated radial velocity 
or proper motion surveys. The future astrometric satellite 
GAIA would be an excellent tool for such a search since it 
will provide accurate proper motions for thousands of stars 
in the halos of nearby globular clusters. 

ACKNOWLEDGEMENTS 
REFERENCES 

Aarseth, S. J., 1999, PASP, 111, 1333 

Allen, C, Moreno, E., & Pichardo, B. 2006, ApJ, 652, 1150 
Baumgardt, H., Makino, J., 2003, MNRAS, 340, 227 
Bekki, K., Couch, W.J., Drinkwater, M.J., 2001, ApJ, 552, 
L105 

Bekki, K., Couch, W.J., Drinkwater, M.J., Shioya, Y., 
2003, MNRAS, 344, 399 

Binney, J., Tremaine, S., 1987, Galactic Dynamics, Prince- 
ton Univ. Press, Princeton, p. 425 

Bruzual, G, Chariot, S., 2003, MNRAS, 344, 1000 

Carraro, G, Lia, C, 2000, A&A, 357, 977 

Casares, J., 2006, in Observational evidence for stellar- 
mass black holes, lAU Symposium 238 in press, 
"arXiv:astro-ph/0612312' 

Capuzzo Dolcetta, R., Di Matteo, P., Miocchi, P., 2005, 
AJ, 129, 1906 

Dabringhausen, J., Hilker, M. & Kroupa, P., 2008, MN- 
RAS, 386, 864 

Dabringhausen, J., Baumgardt, H. & Kroupa, P., 2008, 

MNRAS submitted 
Drinkwater, M.J., Jones, J.B., Gregg, M.D. & Phillipps, S., 

2000, PASA, 17, 227 
Drukier, G, et al., 1998, AJ, 115, 708 

Evstigneeva, E.A., Gregg, M.D., Drinkwater, M.J. & 
Hilker, M., 2007, AJ, 133, 1722 

Fellhauer, M., Kroupa, P., 2002, MNRAS, 330, 642 

Firth, P. et al. 2007, MNRAS, 382, 1342 

Forbes, D. et al. 2008, MNRAS submitted. larXiv :0806. 10901 

Giersz, M., Heggie, D.C., 1994, MNRAS, 268, 257 

Gilmore, G. et al. 2007, ApJ, 663, 948 

Goerdt, T. et al., 2008, MNRAS, 385, 2136 

Hasegan, M. et al., 2005, ApJ, 627, 203 

Heggie, D.C., Hut, P., 2002, The Gravitational Million- 
Body Problem, Cambridge University Press, p. 8 

Hilker, M., Infante, L., Viera, G, Kissler-Patig, M. & 
Richtler, T., 1999, A&AS, 134, 75 

Hilker, M., et al., 2007, A&A, 463, 119 



Jones, J. B et al. 2006, AJ, 131, 312 
Kalirai, J.S. et al., 2008, ApJ, 676, 594 
Kroupa, P., 2001, MNRAS 322, 231 

Kruijssen, J. M. D. 2008, A&A letters in press, 

arXiv:0806.0852 
Lyne, A.G, Lorimer, D.R., 1994, Nature, 369, 127 
Makino, J., Fukushige, T., Koga, M., & Namura, K., 2003, 

PASJ, 55, 1163 
Mashchenko, S., Sills, A. 2005, MNRAS, 619, 258 
McLaughlin, D. E & van der Marel, R. P. 2005, ApJS, 161, 

304 

Mieske, S., et al., 2005, A&A, 430L, 25 
Mieske, S., Hilker, M., Infante, L. & Jordan, A. 2006, AJ, 
131, 2442 

Mieske, S., Hilker, M., Jordan, A., Infante, L. & Kissler- 
Patig, M. 2007, A&A, 472, 111 

Mieske, S., et al., 2008, A&A in press, | aT3Oy:0806.0374l 

Mieske, S., Kroupa, P., 2008, ApJ, 677, 276 

Oh, K.S., Lin, D.N.C., 2000, ApJ, 543, 620 

Pfahl, E., Rappaport, S., Podsiadlowski, P., 2002, ApJ, 573, 
283 

Plummer, H.C., 1911, MNRAS 71, 460 

Rejkuba, M., Dubath, P., Minniti, D., & Meylan, G. 2007, 

A&A, 469, 147 
Scarpa, R., Marconi, G, Gilmozzi, R., & Carraro, G. 2007, 

ESO Messenger 128, 41, arXiv:0707.2459 
Spitzer, L., 1987, Dynamical Evolution of Globular Clusers, 

Princeton University Press 
Thomas, P.A., Drinkwater, M.J. & Evstigneeva, E., 2008, 

MNRAS in press, iirXiv:0801.4840 
Thorsett, S.E., Chakrabarty, D., 1999, ApJ, 512, 288 



8 H. Baumgardt and S. Mieske 



Table 1. Details of the Af-body models. The second column gives the initial number of stars, the third column the initial number of 
dark matter particles. The fourth column gives the mass of a dark matter particle and the fifth column gives the relative mass fraction 
in dark matter and in stars. The last columns give the fraction of dark matter remaining inside the 5% lagrangian radius and inside the 
half-mass radius after one, two and ten apparent relaxation times (eq. IIOI I have passed. 



Nr. 






rriDM 


M, : Mdm 




T = Trh 






JdmIt^ 










[Mq] 






r<RH 


r < R5% 


r<RH 


r < Rsx 


r < R„ 


1 


8224 


91776 


0.030 


1:1 


0.39 


0.78 


0.22 


0.63 


0.01 


0.26 


2 


4285 


95715 


0.030 


1:2 


0.57 


1.51 


0.29 


1.27 


0.03 


0.55 


3 


8223 


27528 


0.100 


1:1 


0.38 


0.81 


0.20 


0.70 


0.01 


0.35 


4 


8224 


183576 


0.015 


1:1 


0.40 


0.78 


0.18 


0.61 


0.01 


0.27 



